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Quadratic  Differential  Equations 

Daniel  E.  Kodltschek  and  Kumpatl  S.  Narendra 
Yale  University 

1.  Introduction; 

Linear  time-invariant  differential  equations  have  been  studied  over  the  past 
three  decades  and  the  stability  properties  of  systems  described  by  the  vector 

e 

differential  equation  x ■ Ax  are  known  to  be  determined  by  the  eigenvalues  of 
Che  matrix  A.  In  contrast,  the  wealth  of  techniques  available  for  nonlinear 
systems  - oscillations,  asymptotic  behavior,  perturbation  theory,  etc.  - reflects 
Che  obvious  fact  that  a broad  generalization  of  one  theory  Is  not  applicable  to 
all  classes  of  systems.  Even  If  our  Inquiries  are  limited  to  the  two  dimensional 
case,  the  volume  of  detailed  work  on  special  systems  such  as  the  Vander  Pol  or 
Dufflngs  equations  and  the  work  of  the  Russian  school  In  the  early  fifties  on 
absolute  stability  demonstrates  the  need  to  treat  these  problems  Individually. 

The  results  of  this  paper  Indicate  that  a profitable  generalization  of  second 
order  differential  equations  may  be  possible  when  we  confine  our  attention  to  Che 
quadratic  case: 

X £ 

and  Its  generalization 


X - Ax  + f(x)  (2) 

which  la  considered  In  a subsequent  paper.  Moreover,  the  Insights  provided  by  this 


initial  study  will  also  aid  In  the  analysis  of  higher  order  quadratic  systems. 


More  chan  being  a convenlenc  class  of  nonlinear  systems,  quadratic  differ- 
ential equations  have  a traditional  Importance  In  stability  literature.  Given 
an  arbitrary  autonomous  dlffsrentlal  squatlon  x ■ g(x),  g(x)  may  be 
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W expanded  in  a Taylor  series  if  it  satisfies  certain  regularity  conditions. 
Lyapunov  proved  that  if  the  equilibrium  state  of  the  linear  approximation  is 
asymptotically  stable  (unstable)  the  nonJ inear  system  will  also  be  asymptotically 
stable  (unstable).  However,  in  the  critical  case  when  the  linear  approximation 
is  merely  stable,  the  higher  order  terms  must  be  examined  to  determine  the  nature 
of  the  stability  of  the  equilibrium  state  and  interest  shifts  to  an  equation  of 
the  form  (2) . 

Recently  the  special  class  of  bilinear  systems 

X • Ax  + u B X (3) 

has  received  a great  deal  of  attention  in  the  control  literature  and  the  princi- 
pal results  of  this  theory  are  influencing  the  direction  of  research  on  general 
nonlinear  systems.  When  the  control  u(t)  in  equation  (3)  is  a linear  function  of 
the  state  variables,  equation  (3)  becomes  a special  case  of  (2)  and  the  stability 
properties  of  such  systems  are  bound  to  be  of  Interest  to  control  theorists. 
Further,  quadratic  differential  equations  are  also  known  to  arise  in  adaptive 
control  [5]  where  the  control  parameters  of  a linear  system  are  continuously  ad- 
justed and  become  state  variables  of  a quadratic  system. 

*'  ^ this  paper  ws  undartal^  a detailed  study  of  the  stability  properties  of 

the  quadratic  differential  equation  (J^  • After  observing  that  such  systenu  can 
never  be  asymptotically  stable  the  equilibrium  states  of  are  classified  in 
terms  of  the  matrices  C and  H.  Necessary  and  sufficient  conditions  for  the  sta- 
bility of  the  origin  are  derived  ia  eeetisw  and  constitute  the  principal  con- 
tribution of  this  paper.  Finally,  these  conditions  are  re-derived  and  elaborated 
using  polar  coordinates  which  allow  a convenient  classification  of  instability 
behavior.  This  exhausts  the  stability  characteristics  (in  the  sense  of  Lyapunov) 
of  second  order  quadratic  differential  equations. 


) 
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2.  Homogeneous  Systems  of  Even  Degree: 


Consider  a dynamical  system  In  R defined  by 


X - h(x)  where  h(cx)  ■ c*Si(x)  and  k > 1. 


The  "homogeneity"  of  h forces  the  direction  of  the  field  to  be  constant  along  any 
straight  line  through  the  origin.  The  consequences  of  this  simple  property  pervade 
the  following  sections. 

If  p(t;xQ)  denotes  the  solution  of  (4)  given  Initial  condition  p(tQ;xQ)  ••  x^, 

k-1 

It  can  be  easily  shown  [2]  that  p(t;xQ)  >•  6p(B  t:xQ).  Now  suppose  k Is  even: 


letting  6 ■ -1  we  have 


p(t;-XQ)  - -p(-t;xQ) 


which  Implies  that  any  trajectory  for  t ^ 0 through  x^  has  an  associated  trajectory 
through  -Xq  for  t ( 0 which  la  Its  reflection.  Hence,  If  n Is  a neighborhood  of 
the  origin  and  ±Xq  ^ d but  +p(T;Xq)e  d for  some  time  T > t^  then  by  (5)  the  trajec- 
tory p(t;-p(T;xQ))  will  leave  0 after  a ^ 

finite  time  and  pass  through  -x_  after  ^ '*"^7 


-p(T;Xo) 


T - t.  time  has  elapsed  [Figure  1].  \ • \ 

° \ p(T;Xo)\ 

This , in  turn,  Implies  that  an  even  homogene-  ^ — \ 

u ^ ^ V,  I -P(T;xo)  I 

ous  system  can  never  be  asymptotically  stable.  I # 

With  this  general  observation  we  set  V J — 

k ■ n ■ 2 and  devote  the  remainder  of  the  Jr 

• 2 ‘Xq 

paper  to  the  study  of  x > f(x)  In  R . 

Figure  1.  Reflection  Property  of  Even 

Instability  of  an  Isolated  Singularity  at  the  Origin:  Homogeneous  System 

The  first  direct  consequence  of  the  homogeneity  property,  (4),  for  the  system 
(1)  Is  that  a second  order  homogeneous  differential  equation  of  second  degree  cannot 
be  stable  unless  its  field  vanishes  on  an  entire  line. 


0 
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Suppose  the  origin  Is  stable:  l.e.  for  any  c > 0 we  can  find  a 6 > 0 such 

that  trajectories  starting  In  B.  - the  ball  of  radius  6 around  the  origin  - always 

c 

remain  In  From  the  preceding  discussion  this  system  cannot  be  attractive. 

Hence,  for  any  p > 0 we  can  find  a v > 0 such  that  trajectories  which  start  outalde 

B never  penetrate  B . Then  by  suitably  choosing  v < ;.  < 5 < e.  any  trajectory 
p V 

starting  In  tl,  • B,  - B will  remain  In  fi_  » B - B for  all  t ^ 0,  and  hence  In  fl, 
lop  zev  ^ 

(see  fig.  2).  By  the  Polncard-Bendlxson 

theory  [3]  any  autonomous  trajectory  In  a 
2 

closed  subset  of  R must  be  either  attracted 
to  a singularity  or  must  be  a closed  path. 

Any  closed  path  trajectory  In  <»u8t 
cross  the  line  through  and  the  origin 
In  at  least  two  different  directions.  But 
this  violates  homogeneity,  hence  ^2  must 
contain  a singular  point.  By  construction, 
this  cannot  be  the  origin.  Since  f(y)  ■■  0, 
y ft  0 Implies  f(ay)  - 0 by  homogeneity,  the 

field  must  vanish  on  the  entire  line  through  this  singular  point  and  the  origin. 

The  next  section  will  Investigate  the  existence  of  such  lines  for  the  quad- 
ratic differential  equation. 

3.  Equilibria  of  Second  Order  Quadratic  Systems; 

We  now  consider  the  particular  class  of  second  degree  second  order  systems 
described  by  (1). 


Figure  2. 


Xj^  ■ X Gx 

• T 
Xj  • X Hx 


(1) 


with  which  this  paper  Is  chiefly  concerned.  We  will  assume  that  neither  G nor  H 
Is  Identically  sero,  and,  without  loss  of  generality,  that  both  are  symmetric. 
From  the  discussion  la  the  previous  section  It  Is  clear  that  the  locus  of  the  set 
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< 


of  critical  points  of  f Is  crucial  to  the  stability  properties  of  (1).  Since 
this  Is  completely  determined  by  G and  H,  well-known  properties  of  symmetric 
matrices  yield  the  following  exhaustive  classification  of  the  types  of  equili- 
brium states  of  the  quadratic  differential  equation.  The  field  In  (1)  may  vanish 

(I)  only  at  the  origin 

(II)  along  a straight  line  through  the  origin 

(III)  along  two  straight  lines  through  the  origin. 

If  either  G or  H Is  definite,  the  system  Is  obviously  of  type  (1).  If  G and  H 

are  Indefinite,  but  of  full  rank,  the  system  may  be  of  type  (1),  (11),  or  (111). 

If  G Is  singular  and  H Is  Indefinite  and  of  full  rank  (or  vice  versa)  the  system 

Is  of  type  (1)  unless  H maps  the  zero  eigenvector  of  G into  Its  orthogonal  comple- 

T 

ment  (l.e.  Gx  ■ 0 & x Hx  « 0)  In  which  case  the  system  Is  of  type  (11).  If  G and 
H are  both  singular  then  the  system  Is  of  type  (1)  unless  G ~ cdl.  In  which  case 
we  get  type  (11).  These  cases  are  Illustrated  In  the  following  examples: 

Example  3.1:  G & H full  rank.  Indefinite. 


(a)  If  G 


' 0 1/2" 

1 0^ 

and  H - 

1/2  0 

0 -1 

, then  X Gx  ■ 


0 Is  satisfied 


by  the  vectors 


and 


T 2 2 

, and  X Hx  ■ Xj  ” *2  “ ® satisfied  by  the  vectors 


and 


1 

-1 


X Gx 


, but]  - I ~ 0 has  only  the  trivial  solution.  The  system  Is  type  (1). 
X Hx 


’ 0 

1/2 

.1 

0 

2 

1 

0 

(b)  If  G - 

1/2 

0 

And  H - G - 

2 

0 

■ 4G  then 

0 

and 

1 

both 


satisfy  X • ^ 


X Gx|  I Xj^Xj 


X Hx  I 1 Ax^Xj 


0 and  the  system  Is  type  (ill). 


€ 


(c)  If  G 


' 0 1/ 

■ 

1/2  0 


and  H • 


1 1/2 


, then  X Gx  » X X ■ 0 for  , 

1/2  0 ^ ^ 0 1 


10 


0-1  0 

, . Hence  x 0 along  and  the  system 

1 1 1 


and  X Hx  - + X2)  “ 0 for  , 

is  type  (11).  L J L 

Example  3.2;  G & H singular. 


1 0 0 0 0 T 2 

(a)  If  G - and  H ■ , then  A satisfies  x^Gx  ” x - 0 

0 0 0 1 1 ^ 


1 T 2 

uniquely,  and  A | satisfies  x Hx  ■ x.  *0  uniquely  so  the  system  Is  type  (1). 

rv  1 i 


10  -10  0 X Gx 

(b)  If  G ■ and  H - -<J  ■■  then  satisfies  „ 

0 0 0 0 1 X Hx 

2->  LJLJ  U 


■ 0 uniquely. 


Example  3.3;  G Indefinite,  full  rank;  H singular. 


(a)  If  G - 


0 1/2 


and  H ■ then  and  satisfy  x*Gx  ■ x,*,  “ 0 

1/2  0 11  0 1 ^ ^ 


1 1 


T 2 2 

and  satisfies  x Hx  ■ x^^  + ^*1*2  ^ *2  " ® uniquely,  hence  the  origin  Is  the 


only  equilibrium. 


loriuB.  — 

0 1/2 

If  G - 

1/2  0 


0 0 1 

and  H ■ then  , the  zero  eigenvector  of  H, 

0-1  0 


Is  mapped  by  G Into  - Its  orthogonal  complement.  The  system  Is  of  type  (11). 

[1/2  J 

In  sunnary,  the  necessary  condition  of  a zero  line  stipulates  that  stable 
systems  must  be  of  the  kind  3.1.b,  3.1.C,  3.2.b,  3.3.b.  Sufficient  conditions 


are  derived  In  the  next  section. 
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4.  Stability  of  Second  Order  Quadratic  Systems; 


From  the  discussion  in  section  2 it  follows  directly  that  any  quadratic  sys- 
tem of  type  (1),  with  a unique  isolated  equilibrium  at  the  origin,  is  unstable. 

In  this  section  we  shall  further  refine  necessary  conditions  and  finally  arrive 
at  conditions  which  are  both  necessary  and  sufficient  for  stability.  This  is  most 
easily  accomplished  by  first  deriving  a special  form  for  f(x)  in  (1)  that  charac- 
terizes systems  of  type  (11)  and  (ill). 

Quadratic  Differential  Equations  with  Non-isolated  Equilibria; 

If  the  origin  is  not  an  Isolated  equilibrium  (l.e.  if  the  system  is  not  of 

T T 

type  (1))  then  let  lleil  - 1 and  f(e)  - 0 (or  e Ge  » e He  *=  0) . By  homogeneity  (4) 
we  know  the  whole  line  Xe  (X  e IR)  is  a set  of  critical  points  of  f.  Now  define  a 
new  coordinate  system  by  the  orthogonal  transformation  y = R x where  R = [e.c] 


(Kell  ■ l|c||  ■ 1 and  e c ■ 0) . Then 


where 


T T 

c Ge  c Gc 


f(x)  - f(Ry)  - 


and  Q - 


T 

0 e He 

T T 

c He  c He 


Hence  we  can  factor  a y2  ” c x out  of  both  quadratic  forms  to  write 


where  D ■ 


r T T ' 

2e  Gc  c Gc 

■ j .p 

2e  He  c He 


f(x)  - c x Dx 


c*G  2 0 


0 1 


Equation  (7)  expresses  the  fact  that  any  quadratic  system  whose  field  vanishes 


along  at  least  one  line  has  solutions  which  are  associated  with  a linear  system 


z - BDz  (0  e R) 


I 

i 

( 
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Glven  a fixed  initial  condition,  z , the  solution  of  (7i)  for  any  constant  0 is 

° 0 

z„(t)  ■ , thus  Zg  (t)  “ Zg  (t)  where  t - t.  In  other  words,  the  trajec- 

2 

torles  of  (7i)  cut  out  identical  manifolds  in  R for  each  0,  although  their  time 
parametrlzation  on  these  manifolds  is  dilated  or  contracted  by  an  appropriate 
constant.  More  generally,  if  we  are  given 


then  x(t)  “ e‘ 


0(T)dT 


0* 


X - 0(t)Dx  (7t) 

Again,  the  solution  of  (7t)  for  any  scalar  time  function 


0(t)  given  initial  condition  x will  be  contained  in  the  same  manifold  as  z (t); 


however,  its  time  parametrlzation  will  be  varied  by  the  function 


0(T)dT,  Thus,  If 


T 

we  think  of  (7)  as  a particular  case  of  (7t)  where  0(t)  » c p(t;xQ)  has  been  com- 
puted a priori,  then  we  see  that  all  trajectories  of  the  quadratic  differential 
equation  with  a non-lsolated  set  of  equilibria  must  lie  on  manifolds  determined 
by  Che  linear  system  (71).  With  this  observation  in  mind,  we  may  now  extend  neces- 
sary conditions  for  stability  by  considering  different  classes  of  linear  second 
order  time-invariant  systems.  Specifically  we  will  examine  the  nature  of  (7)  when 
the  equilibrium  state  of  the  linear  system  (71)  is  a node,  center  or  focus, 
a.  D singular; 

T 

Suppose  D ■ kb  - a singular  matrix.  Then  (7)  may  be  re-written  as 


X - |^x’'bc^xj  k (8) 

T 

The  symmetric  part  of  be  is  indefinite  when  b c.  In  this  case  (8)  is  of  the 

* T 

form  X ■ k X Gx,  hence  G ■ **  indefinite  and  from  the  discussion  in  section 

3 this  is  the  only  possible  occurrence  of  a system  type  (ill).  Clearly  this  equa- 
tion defines  a field  whose  direction  is  uniformly  specified  by  the  vector  k and 
which  vanishes  on  the  two  lines  Xe^,  Xe2  (Xc  IR) , satisfying  e^  Gk  0 
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then  either  p(t;k)  or  p(t;-k)  will  tend  to 

Infinity  along  the  vector  k (or  -k)  as 
T 

t o since  X Gx  has  the  same 

X - X(t)k 

sign  for  t i t^.  This  is  illustrated  in 
figure  3. a. 

A more  interesting'  case  arises  when 
T 

x Gx  0 is  satisfied  on  the  line  determined 

by  k.  Let  r also  satisfy  this  equality 

(r  ^ k)  and  choose  the  four  points  Figure  3. a.  Type  (ill):  two  zero  lines 

unaligned  with  field  direction. 

p - Ikiar  (a  e R) . The  field  vector  at  these  points  is  either  in  the  positive  or 

T T T 

negative  direction  of  k as  determined  by  p Gp  » 2ak  Gr  or  -2ak  Gr.  If  a > 0 and 
T T 

we  assume  k Gr  < 0 (note  k Gb  0 since 

k ft  b)  it  follows  immediately  that 

p(t;k  - or)  “ in  the  direction  of  k, 

and  p(t;-k  - or)  -»  “>  in  the  direction 

of  negative  k as  shown  in  figure  3.b. 

T 

If  k Gr  < 0 then  the  alternate  pair  of 
trajectories  is  unstable.  This  exhausts 
the  possible  occurrences  of  type  (ill) 
equations.  Hence  necessary  conditions 
for  stability  have  been  extended  to  ex- 
clude all  but  systems  of  type  (il)  - with 

Figure  3.b.  Type  (ill):  field  direction 
a single  line  of  equilibria.  aligned  with  one  of  two  zero  lines. 


The  remaining  case  of  singular  D occurs  when  the  matrix  be  in  (8)  is  symmetric: 

i.e.  b ■ c.  If  k c^  then  either  p(t:k)  or  p(t;-k)  follows  the  line  determined  by 

T T 

k to  infinity.  If  k ■ c^  then  equation  (7)  may  be  written  x - c^  x cc  x and  every 
trajectory  which  has  non-zero  motion  tends  to  infinity  in  the  direction  c^  as  shown 
in  figure  4. 
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b.  D Non-Singular  with  Real  Eigenvalues: 

If  the  equilibrium  of  (71)  is  a node 
(i.e.  Che  eigenvalues  of  D are  real  and 
non-zero)  then  its  associated  quadratic 
system  (7)  must  be  unstable.  If  D has 
an  eigenvector,  k,  which  is  not  orthogonal 
to  c then  either  p(t;k)  or  p(t;-k)  must 
tend  Co  infinity  in  the  direction  k.  This 
is  illustrated  in  Figure  5. a.  If  no  such 


Figure  4.  Type  (ii) : field  direction 
aligned  with  unique  zero  line 


k exists  then  D must  be  non-dlagonalizable 


Figure  5. a.  Type  (11):  Pure  Nodal 

Behavior 


and  c^  its  only  direct  eigenvector.  In 
this  case,  since  all  manifolds  of  (7t) 
Intersect  the  equilibrium  set  of  (7)  only 
once  (excluding,  of  course,  that  set  itself) 
as  depicted  in  Figure  5.b,  and  since  the 
solutions  on  these  manifolds  are  mirror 
reflections  across  the  origin  (from  section 
2),  trajectories  on  one  side  of  the  zero 
line  of  (7)  must  tend  to  infinity. 

By  this  argument  we  have  excluded  all 
systems  of  form  (7)  whose  D matrix  has  real 
eigenvalues.  Thus  the  following  is  a state- 
ment of  the  refined  necessary  conditions: 
the  quadratic  differential  equation  (1)  is 
stable  only  if  its  field  vanishes  on  a 
single  line  through  the  origin  and  its 
associated  linear  system  (71)  exhibits 


Behavior 

behavior.  We  will  now  demonstrate  Chat  this  condition  is  also  sufficient. 
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c . D with  Complex  Conjugate  Eigenvalues: 

If  the  solution  manifolds  of  (71)  Intersect  every  direction  on  the  plane  at 
a finite  point  at  least  once  then,  recalling  the  condition  from  section  2 that 
every  trajectory  of  (1)  must  be  contained  in  a half  plane  we  might  suspect  that 
(7)  is  stable  when  D has  complex  conjugate  eigenvalues.  In  this  case,  the  con- 
taining half  plane  Is  defined  by  the  zero  line  c^  and  the  solution  to  (7)  is 

trapped  on  a half  loop  of  the  spiral  or  circle  defined  by  the  solution  of  (71). 

This  can  be  argued  more  formally.  Define  ^(xq)  “ e (t  < “) (y  ■ e^^x^)) 

- the  manifold  cut  out  by  a particular  solution  of  (71)  given  initial  state  x^, 
and  the  half  plane  P = {y  e R ly  ^ 0^*  If  pCiJXq)  Is  a solution  of  (7)  then 
(^t)  p(t;xQ)  e L(xQ)nP.  Since  (1)  is  well  behaved,  the  fact  that  L(xQ)nP  may  be  dis- 
connected is  of  no  concern-we  assume  that  p(t;xQ)  lies  entirely  on  the  branch.  A,  of 

PnL(xQ)  which  is  connected  to  Xq.  Clearly  A is  contained  within  the  closed  half-annulus 
of  radius  HxqII  ±6  for  some  constant  6.  Having  ruled  out  the  possibility  of  closed 

f\j  . 

paths  we  must  have  (by  Bendixson  Theory)  p(t;xQ)— >x,  x e PnLCx^).  Inspecting  (7) 
it  is  obvious  that  x lies  on  the  line  orthogonal  to  c.  Hence  for  an  arbitrarily 
small  e > 0 we  can  always  choose  |Ixq||  small  enough  to  obtain  p(t;xQ)  e B^(^t  > Iq)  • 

The  system  is  stable.  Since  (1)  is  homogeneous  any  stability  characteristic  must 
be  global.  Hence  we  have  the  following: 

Theorem:  The  second  order  quadratic  differential  equation  (1)  is  globally  stable 
if  and  only  if  it  is  of  the  form 

. T 

X ■ c xDx 


Example;  If  G ■ 


and  D 


les  of  thcL  matrix  D 

are  complex. 

p -j 

0 1/2 

-.76  -.is" 

. T 

1 

and  H * 

then  X ■ c xDx  with  c ■ 

1/2  0 

-.15  0 

0 

L J 

-.76  -.30 


Since  Che  eigenvalues  of  D are  complex  Che  system  is  stable. 


The  system  was  simulated  and  results  are  plotted  in  Figure  6. a. 

0 1/2 

-1  0 

• T 

1 

Example:  If  G ■> 

and  H ■ 

then  X ■ c xDx  with  c ■ 

1/2  0 

0 0 

.0- 

i 


Figure  6a.  Analog  Computer  Simulation  of  a Stable 
Quadratic  System  with  Focus  Behavior: 
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and  D ~ 


Since  the  eigenvalues  of  D are  Imaginary  the  system  Is  stable. 


0 1 
L-i  0 

Tha  system  was  simulated  nnd  results  are  plotted  In  figure  6.b. 

5.  Global  Behavior  of  Second  Order  Quadratic  Systems; 

Results  of  earlier  sections  Indicate  that  most  quadratic  systems  are  unstable. 
Hence  the  derivation  of  necessary  and  sufficient  conditions  for  stability  cannot 
address  the  qualitative  behavior  of  the  solutions  In  most  Instances  of  equation  (1). 
To  fill  this  gap,  and  since  the  Investigation  of  Instability  behavior  Is  an  Intrin- 
sically Important  component  of  any  stability  analysis,  this  section  will  be  con- 
cerned with  the  classification  of  the  global  properties  of  second  order  quadratic 
systems. 

a.  Polar  Coordinate  Representation: 

For  the  purposes  of  this  Investigation  It  Is  most  convenient  to  express  equa- 

A n 2 ^^2 

tlon  (1)  in  polar  form.  Using  the  transformation  p •=  /x  + x„,  0 ■ arc  tan  — we 

X 

may  write 

p - 1/p  (x^Xj^+X2*2^ 

a 2 * * 

0 •»  1/p  (x^X2-X2Xj^) 

• T • T A 2 

where  x^  - x Gx  and  X2  - x Hx  from  (1).  If  we  define  g(v)  - g^v  + g^v  + g^^  and 

A 2 ^ 

h(v)  - h^v  + hpV  + hj^,  and  let  0 e (-it/2,ir/2)  then  this  may  be  re-wrltten  as 

• 2 3 

p - p cos  0[g(tan  0)  + tan  0 h(tan  0)] 

3 

0 ■ p cos  0[h(tan  0)  - tan  0 g(tan  0)) 

Recalling  from  section  2 that  the  solutions  of  (1)  behave  symmetrically  across  any 
line  through  the  origin.  It  la  clear  that  there  Is  no  loss  of  generality  in  re- 
stricting 0 to  the  giver,  open  Interval. 


Let  G - 


1/2  gQ 

and  H - 

■ •'l 

1/2 

1/2  gj 

*2 

[l/2  h„ 

*'2  _ 

1 
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Once  again,  the  acablllty  analysis  of  (9)  requires  an  Identification  of  Its 

• • 

equilibrium  states.  If  p = 0,  then  p ■>  0 - 0;  the  abscissa  In  the  new  coordinate 

system  (the  origin  in  (1))  Is  an  equilibrium  set.  For  p ¥ 0,  the  equilibrium 

states  of  (9)  are  completely  determined  by  0,  which  Is  a direct  consequence  of 

A *2 

homogeneity.  Parametrizing  the  open  right  half  plane  by  the  scalar  v ■>  tan  © " — , 

*1 

the  critical  points  of  (9)  are  obtained  by  simultaneously  solving  the  equations 

r(v)  ^ g(v)  + vh(v)  ■ 0 
p(v)  ■ h(v)  - vg(v)  » 0 

where  r or  p (or  both)  are  cubic  polynomials  in  v.  In  the  sequel  we  will  assume 
that  p Is  cubic:  l.e.  g2  0 with  no  loss  of  generality.  (If  g2  ■ 0,  a new 
Cartesian  system  may  be  chosen  by  an  orthogonal  transformation.  This  Is  equiva- 
lent to  the  choice  of  a different  parametrlzatlon  of  0).  Thus  p(v)  ” 0 must  have 
either  three  real  roots,  two  real  roots  (one  with  multiplicity  2) , or  a unique  real 


root. 

The  following  result  summarizes  the  character  of  solutions  of  (10),  and,  hence, 
the  equilibria  of  (9). 

Lemma:  Vq  Is  a real  root  of  (10)  If  and  only  if  It  Is  a real  root  of  the  simultane- 
ous equations  g(v)  ■ 0,  h(v)  ■ 0. 

Proof : The  forward  direction  Is  trivial.  To  prove  the  converse,  asstime  that  v^  Is 
a real  root  of  (10)  (l.e.  rCv^)  - 0 • p(Vq)),  but  g(vQ)  0.  In  this  case  we  have 
-1 

—7 — - — : l.e.  V-  must  be  purely  Imaginary  - contradiction. 

g(vQ)  0 ° V , t.  j 

By  this  lemma,  (10)  must  have  either  (1)  no  solutions,  (11)  one  solution,  or 
(111)  two  solutions,  with  the  following  characteristics: 

(1)  If  equation  (10)  has  no  solutions,  system  (9)  has  the  trivial  equilibrium 


set  p ■ 0.  However,  since  p is  cubic  there  must  be  at  least  one  real  v^  such  that 
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h(vQ) 

p(v.)  - 0.  For  any  such  roots  the  relation  —7 — r “ v.  must  be  true. 

0 ^ 8(Vq)  0 

(II)  If  equation  (10)  has  a single  solution  Vq,  then  system  (9)  has  an  equili- 

brium set  along  the  single  line  0 - arc  tan  v^.  If  p has  other  distinct  roots,  then 
they  cannot  be  roots  of  g by  the  lemma:  hence  the  relation  ^ must  be  satis- 
fied at  these  points.  If  Is  the  unique  real  root  of  p then  the  relation  ^ 

Is  never  satisfied  for  v c R.  (Unless  has  multiplicity  3 in  p). 

(III)  If  equation  (10)  has  two  distinct  solutions  and  then  system  (9) 

has  equilibria  along  two  lines  0 *•  arc  tan  v and  6 ■ arc  tan  v*.  Since  p Is  cubic 

» » h(v;) 

and  must  have  a third  root,  vJJ,  which  cannot  solve  (10),  the  relation 


will  hold. 

It  Is  clear  that  these  three  cases  correspond,  respectively,  to  the  three  types 
of  equation  (1)  Identified  In  section  3.  Using  the  results  of  section  4 we  may 
Immediately  see  that  the  necessary  and  sufficient  condition  for  stability  is  given 
by  the  second  alternative  of  case  (11).  This  may  be  confirmed  Independently  as 
follows. 

Stability  Analysis  of  Polar  Form: 

Assume  case  (11)  holds  and  v^  Is  the  unique  real  root  of  p(v)  ■ 0.  Then  for 
some  0 < M < «,  -M  $ p(^  ^ fro®  (9),  we  have  p « Q i pM0.  Since 

p ^ 0 we  may  write 


(11) 


Now,  factoring  p(v)  ■ (v-Vq)p'(v),  we  must  have  p'(v)  bounded  away  from  zero. 

If  p'  Is  always  negative  let  ^ ^ 0 - 0^,  If  p'  Is  positive,  then  let  ^ - 0 - 0^  - w. 
Then  (assuming  the  former)  we  have 


p(t)  co8^(4+0Q)p'Itan((>46Q)]j(tan(^+0Q)  - v^)  cos(#+0q^ 

or,  multiplying  the  polynomials  In  tan  0 through  by  cos  0 

- (sln(ffeQ)  - Vp  cos (446^)1  a (t,0)  (12) 

♦e  (-0Q-ir , -Oq+v) 


I 
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where  a(t,0)  Is  a negative  definite  function.  Since  (>[8ln((K6Q)  - co8(^+©q)]  i 0 

equation  (12)  Is  an  asymptotically  stable  scalar  system,  and  4i(t)  -*■  0 for  all  Initial 
conditions  ♦p  e (-0^-*  .-O^+Tr) . 

Since  In  p la  a positive  function,  the  stability  of  (12),  and  the  majorlzatlon 

In  (11)  forces  the  boundedness  of  In  p.  Of  course  In  p la  bounded  If  and  only  If  p 

la  bounded.  Then  If  - v does  not  hold  for  any  v,  equation  (1)  Is  globally 

g(v) 

stable. 

Conversely,  It  may  be  similarly  demonstrated.  Independent  of  the  analysis  In 
section  4 that  all  remaining  cases  of  equation  (9)  are  unstable.  Recall  that  In 


each  remaining  case  there  exists  at  least  one  real  v.  such  that 


g(VQ) 


v^  and 
0 


r(vQ)  ^ 0.  Then  for  Initial  conditions  (p»®q)  when  0^  • arc  tan  v^,  system  (9) 


may  be  re-vrltten  as: 


• 2 3 

p - P cos  0Q  r (tan  0^) 
0-0 


which  Is  a scalar  quadratic  system.  Since  p ^ 0,  equation  (8)  has  bounded  solutions 

3 3 

If  and  only  If  cos'*  0^  r (tan  0^)  < 0.  However  If  this  Is  the  case  then  cos  (0^  + w) 

r[tan(0-  + n)l  > 0 and  p(tan  (0-  + x))  - 0.  Hence  Initial  condition  (p  0.  + ir) 

0 0 » 

also  yields  equation  (13)  and  has  the  solution 


* -(2  co8(0q  + n)tan  0^)1  - p^ 

which  escapes  to  Infinity  In  finite  time. 

In  susnary,  we  have  re-stated  the  necessary  and  sufficient  conditions  for  sta- 
bility of  equation  (1)  In  polar  form:  system  (9)  Is  stable  If  and  only  If  there  Is 

h(Vo) 

no  Vq  e IR  such  that  - v^. 


The  crucial  equation  ■■  v Is  Intimately  related  to  the  quadratic  operator 

In  (1) : the  existence  of  a real  to  solve  that  equality  Is  equivalent  to  the 
* Excluding  the  case  where  the  equality  holds  at  v^.and  v^  Is  simultaneously  a 
triple  root  of  p. 
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I 


I 


l 


existence  of  a non-zero  fixed  direction  of  f. 


Suppose  0,  Xq  0).  We  have 


But 


T 

X Mx 
X Gx 


’*2 

h(-^) 

X 


Xx 


02 


Xx 


A ^ 

i-  . Hence  letting  v ■ — ^ we  have 


*2 

g(;^) 

*1 


0 X 


01 


rCVq) 


01 


''o- 


02 

‘01 


The  converse,  while  strongly  suggested  by  the  condition  p 0,  0 ■ 0 is  inmediate 
if  we  recall  the  identity 


' x(0)  - Xq 

2 

where  X2(Xj^)  is  a local  parametrlzatlon  of  the  phase  curve  {p(t;xQ)|t  > 0}  in  IR  . 

Naaely,  the  ratio  of  the  time  derivatives  specifies  the  slope  of  the  tangent  line  of 

h(v  ) 

— ; — T •>  V-  is  equivalent 
rCVq)  0 

on  the  line  X2  - 
, we  must  have  X 0. 

Necessary  and  Sufficient  Conditions  for  Stability; 

In  sunnary  we  have  shown  that  the  following  statements  are  equivalent 

(i)  Equation  (1)  is  globally  stable. 

T 

(11)  f(x)  ••  c xDx  and  D has  pure  complex  eigenvalues. 


the  manifold  on  which  a solution  lies.  Then  the  condition 

I lies 


dx2  X 

to  the  condition  3 — ■ — , and  trajectory  p[t; 
ox . X _ 

v-t- 


Hence  f 


.1 

a 

av 


Oj 


and  since  ^ 0 implies  p ^ 0 


Stability; 


(ill)  The  equation  ^ has  no  real  solutions. 

(iv)  The  quadratic  operator  f has  no  non-zero  fixed  directions. 


and  are  mutually  exclusive  of  the  set  of  equivalent  statements 

** 

(1)  Equation  (1)  has  a finite  escape  trajectory. 


(il)  f is  type  (i)  or  D has  real  eigenvalues. 

Instability; 

(ill)  There  exists  a v.  c R such  that  —7 — r • 

0 g(v-) 


** 


(iv)  The  quadratic  operator  f has  a non-zero  fixed  direction. 

* Excluding  the  case  where  the  equality  holds  at  and  v^  is  simultaneously  a 
triple  root  of  p.  (see  next  footnote)^ 

**  Excluding  the  situation  where  f(x)  ■ c xDx,  D has  a double  real  eigenvalue,  and 
the  only  eigenvector  of  D Is  on  the  orthogonal  complement  of  c.  This  Is  unstable, 
as  shown  In  section  4.  It  Is  equivalent  to  the  condition  that  Vq  Is  a triple  root 
of  p In  case  (11). 


I 
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b.  Classification  of  Unstable  Quadratic  Systems; 

Hsvlnii  rslterated  earlier  results  using  the  polar  coordinate  transformation, 
we  will  now  classify  the  various  types  of  Instability  behavior  solutions  Co  syaton 
(1)  awy  evince  on  the  basis  of  their  associated  phase  curves. 

It  has  already  been  established  that  the  solution  manifolds  of  type  (11)  and 

T 

(111)  systems  given  by  equation  (7)  (l.e.  where  f (x)  ■ c xDx)  are  Identical  to  those 
which  are  associated  with  the  linear  system  (71). 

z • BDz  (71) 

It  Is  well-known  that  system  (71)  may  possess  nodal,  focus,  or  center  characteris- 
tics. In  the  case  of  focus  or  center  behavior  of  (7l)  we  have  demonstrated  that 
(1)  Is  stable.  Otherwise  the  trajectories  of  (7)  follow  half  the  nodal  manifolds 
(71)  to  the  origin,  and  the  other  half  out  to  Infinity,  with  trajectories  starting 
on  an  eigenvector  of  D escaping  to  Infinity  In  finite  time,  (Figures  3-6]. 

On  the  other  hand,  we  have  said  very  little  about  the  trajectories  of  systems 
of  type  (1).  Again,  It  la  more  convenient  to  examine  the  solutions  of  (9)  than  to 
attempt  an  Immediate  analysis  of  the  Cartesian  system  (1).  The  tangent  slope  of 
phase  curves  In  the  (p,0)  plane  Is  given  by 


q(v) 


m ^ 

” de 


r(tan  9)  _ g(v)  l-  vh(v) 

^ p(tan  0)  ” h(v)  - vg(v) 


(12) 


The  poles  of  q correspond  to  the  values  of  v such  that  » v as  established  before. 

The  significant  characterlatics  of  Che  manifolds  of  (1)  are  determined  quite  easily 
by  a glance  at  the  sign  of  q in  the  neighborhood  of  a pole.  If  is  a single  root 
of  p(v)  ■ 0 then  q(vQ  c)q(vQ  - c)  < 0.  If  qfv^  + e)  > 0 then  solution  manifolds 
on  Che  (p,0)  form  a negative  cusp  around  the  line  0 ■ arc  tan  as  shoim  in 
figure  7. a.  If  qCv^  -f  c)  < 0 Chen  solution  oumifolds  on  the  (p,6)  plane  form  a 
positive  cusp  around  the  line  e • arc  tan  v^,  as  shoim  in  figure  7.b.  If  is  a 


double  root  of  p(v)  ■ 0 than  4-  c)q(vQ  - e)  > 0 and  the  solution  manifolds  will 


j 


have  positive  (or  negative)  slope  on  both 
sides  of  the  line  6 - arc  tan  v^  as  shown 
In  figure  7.c. 

Loops  and  Flows; 

In  the  light  of  the  above  observations, 
a rough  translation  of  figures  8 Into  the 
Cartesian  plane  convinces  us  that  the  proto- 
typlc  behavior  of  system  (1)  type  (1)  falls 
Into  two  distinct  classes.  The  trajectories 
must  be  separated  by  at  least  1 and  up  to  3 
separatrlces: the  fixed  directions  of  the  map 
f as  determined  by  the  zeros  of  p.  If  f has 
a unique  fixed  direction,  then  the  solutions 
of  (1)  must  correspond  to  either  figure  7. a 
or  7.b.  In  the  (Xj^.x^)  plane  these  proto- 
types are  given  by  figure  8. a or  8.b  which 
we  will  call  the  "loop"  or  "flow"  prototypes 
rcspecClTely.  If  f has  two  fixed  directions, 
then  p must  have  one  double  root.  In  this 
case,  one  fixed  direction  will  have  either 
the  loop,  or  the  flow  characteristic,  and 
the  other  will  heve  the  special  feature 
depleted  In  7.c.  The  solutions  In  the  (x^,X2) 
plane  retain  either  a flow  or  a loop  character 
as  shown  In  figure  8.c  (flow  dominant)  and 
f.d  (loop  dominant).  Finally,  If  f has  three 
fixed  directions,  p must  have  three  roots  and 
each  root  must  evince  either  pure  loop  or  flow 


q|Vj,-€)>0 

q{v_+€)>0 


q (Vo-€)<0 


(b)  NEGATIVE  CUSP 

q (vo-<)  >0 


(c)  MIXED 

Figure  7.  Characteristic  Manifolds 
on  the  Polar  Plane 


charac  ter la  t Ice 


Illustrative  examples  of  these  cases 
are  now  given: 

Example  5.1. a.  Flow  with  one  fixed  direction: 


Let  G - 


and  H 


g(v)  ■ V ,h(v)  - 1 and  r(v)  ■ v(v  + 1), 

3 

p(v)  ■ 1 - V . The  system  Is  tjrpe  (1) 
since  g and  h have  no  common  factors.  The 
cubic  polynomial  p has  a single  root  at  v • 1, 

r^i 

hence  the  direction  Is  Invariant  under 


f.  The  ratio  g(v) 


a)  Flow 


Is  positive  In 


the  left-hand  neighborhood  of  1.  Hence  the 
tangent  lines  ^ are  given  by  figure  l.b. 

The  actual  system  was  simulated  on  the  digital 
computer  and  results  are  plotted  In  figure  9. a. 
Example  3. l.b.  Loop  with  one  fixed  direction: 


b)  Loop 


Let  G ■ 


1/2  -1 


and  H ■ 


-1  1/2 


2 2 

Then  g(v)  - -(v  -v-1) , h(v)  • v +v-l)  and 
3 3 

r(v)  - V +1,  p(v)  » V - 1.  Again  the 
system  Is  type  (1) , and  p has  a unique  root 


V ■ 1.  The  vector 


Is  again  fixed  under 


f.  Now,  however,  the  ratio  q(v)  Is  negative 

dp 

In  the  left-hand  neighborhood  of  1.  Hence  ^ 

Is  given  by  l.a.  The  actual  system  was  simulated 
on  the  digital  computer  and  results  are  plotted 
in  figure  9.b. 


c)  Mixed  - (Flow  Dominant) 


d)  Mixed  - (Loop  Dominant) 

Figure  8.  Flow  and  Loop  Characteristics 
of  Trajectories 


Unstable  Quadratic  System  with  Loop 
Behavior : 


t 
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c.  The  Stability  of  Off-Orlgln  Equilibria: 

In  Che  preceding  discussion  a polar  coordinate  transf onnatlon  elucidated  the 
analysis  of  the  global  behavior  of  equation  (1)  with  respect  to  the  origin.  To 
complete  the  Investigation  of  global  behavior  an  account  of  the  stability  properties 
of  other  equilibria  are  presented  here. 

Systems  of  type  (1) , with  no  other  equilibria  than  the  origin  have  been  treated 
In  the  previous  subsection.  For  the  two  remaining  types,  recall  from  section  A 
that  equation  (1)  may  be  expressed  using  equation  (7) 

• T ^ 

where  D Is  either  of  full  rank  or  singular.  In  the  former  case  the  system  has  a 
unique  line  of  equilibrium  i.e.  x “ c^^.  In  the  latter  case  the  system  has  two  lines 
of  equilibria  which  satisfy  x =•  c^  and  Dx  = 0 [a  special  case  occurs  when  these  two 
lines  coincide]. 

Systems  of  Type  (111): 

If  D Is  singular,  It  was  shown  In  section  4. a that  equation  (7)  can  be  rewritten 
as  equation  (8): 

• T 

X ■ kx  Gx 


This  is  type  (111)  if  and  only  if  G Is  indefinite  and  of  full  rank.  As  explained 
in  that  discussion  solutions  lie  on  parallel  lines  in  the  direction  k.  These,  in 
turn,  correspond  to  the  solution  manifolds  of  a singular  linear  system  whose  unique 
non-zero  eigenvector  Is  k.  Superimposed  on  this  field  Is  the  equilibrium  set  speci- 


fied by  the  two  distinct  vectors  e^^  and  e^  for  which  x Gx 


- 0 (1  « 1,2)  is 


[x-ei 

satisfied.  In  appendix  1 we  present  a formal  analysis  of  the  perturbed  motion 
around  the  equilibria  Xe^  and  Ae^t  A c 

If  k ^ e^  (1  " 1,2)  then  either  positive  or  negative  open  half  line  (I.e.  ex- 
cluding the  origin)  contains  stable  equilibria  and  its  opposite  closed  half  line 
(i.e.  including  the  origin)  contains  unstable  equilibria.  Each  stable  half  line 


Is  a locally  attractive  Invariant  set.  This  situation  is  portrayed  in  figure  3. a 
(section  A) . 

If  k ■ e^^  then  the  direction  of  the  field  Is  Itself  a line  of  equilibria,  all 
of  which  must  be  unstable.  In  this  case,  the  remaining  line  of  equilibria  has  an 
open  half  line  of  stable  and  a closed  half  line  of  unstable  points.  The  unique 
stable  half  line  is  a locally  attractive  invariant  set.  This  situation  is  portrayed 
in  figure  3.b  (section  A). 

Systems  of  Type  (ii) ; 

As  established  in  section  3 if  D is  nonsingular  (or  singular  and  possesses  a 
zero  eigenvector  in  the  direction  orthogonal  to  c)  then  the  system  is  type  (ii). 
Again  the  solutions  lie  on  manifolds  determined  by  the  associated  linear  system 

z •»  Dz 

Interrupted  by  the  zero  line.  In  appendix  1 it  is  formally  shown  that  in  all  cases 
there  must  be  a half  line  of  stable  equilibria  (excluding  again  the  case  where  the 
only  eigenvector  of  D is  orthogonal  to  c) . The  half  line  is  closed  if  and  only  if 
D has  purely  complex  eigenvalues.  This  line  is  a locally  attractive  invariant  set. 
Two  representative  systems  are  depicted  in  figures  5 


(section  A). 


« 


t 
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Appendlx  I;  The  Stability  of  Off-Origin  Equilibria 

Recall  that  system  (1),  x • f(x),  may  be  re-written  as  equation  (7), 

• T 

x ■ c xDx,  if  there  are  any  off-origin  equilibria.  Around  an  arbitrary  equili- 
T 

brium  state,  Ae(e  c - 0 and  X e |R) , the  perturbed  equation  of  motion  z ■ f(e  + z) 
- f(e)  may  be  simplified  and  written  as: 

• T T 

z • c zDXe  + c zDz 


Define  the  matrix  R - (Dele]  and  the  change  of  basis  y ■ R z. 

• *1.  T T T T 

Then  y • XR  Dec  Ry  + R DRyc  Ry.  Since  c R = [c  De,0]  and  R De 

re-write  the  last  equation  as 

■ ^ T 

y ” De 


I we  may 

oj 


(A.l) 


T 

Fact : System  (A.l)  is  locally  stable  if  and  only  if  Xc  De  < 0 in  which  case  the 

line  y^^  ■ 0 is  an  attractive  set. 

T 2 • T 2 2 

Proof : If  Xc  De  > 0 let  v • y^^.  Then  v =*  Xc  Dey^^  + y^^OCy  ).  Since  v > 0 in  the 

region  v > 0 for  small  enough  values  of  lly|| , the  system  is  unstable  by  Chetaev’s 

theorem. 


T T • T 2 2 

If  Xc  De  < 0,  let  v » l/2y  y.  Then  v * Xc  Dey^^  + y2^0(y  ) lAlch  is  negative 

semi-definite  in  a small  enough  neighborhood  of  the  origin, 0.  Since  v la  always 

on  {Xj^  ^ 0}nn  we  know  {x^^  - 0}nfJ  is  attractive.  Since  the  equation  of  perturbed 

motion  about  any  point  in  {x^  > O)  is  of  the  form  (A.l),  the  same  argument  works 

on  the  entire  line. 

If  the  system  is  type  (ill)  then  using  equation  (8)  we  know 


X - kx^Gx  - kx^bd^x  (8) 

Thus  the  matrix  D in  (7)  may  be  either  kb^  or  kd^  with  c “ d or  c ■ b respectively. 

T 

If  k j*  d^  and  k b^  then  Xc  De  < 0 for  either  X < 0 or  X > 0 using  either  form  of 

Of  course  this  Is  not  possible  if  De  ■ Xe.  In  this  case,  either  choose  the 
second  zero  line  (If  the  system  is  type  (ill))  and  perform  the  Identical 
analysis,  or  If  the  system  is  type  (11)  then  there  Is  no  attractive  half  line, 
(see  second  footnote  on  p.  16). 
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« 


D.  Thus  chere  are  two  attractive  open  half  lines.  If  k^-  d or  k^-  b then 
T 

c De  ■ 0 for  one  choice  of  D.  Hence  there  Is  only  one  open  attractive  half  line. 

If  the  system  Is  type  (11)  then  the  matrix  D In  equation  (7)  Is  unique  and 

T T 

either  nonsingular  or  equal  to  cc  (see  section  4).  If  c De  0 then  there  must 

be  an  attractive  half  line  (open  or  closed  depending  upon  the  stability  of  the 

T * 

origin  of  course).  Notice  c De  ■ 0 only  If  e is  the  unique  eigenvector  of  D 

(see  footnote  on  previous  page). 


Of  course  this  Is  not  possible  If  De  - Ae.  In  this  case,  either  choose  the 
second  zero  line  (If  the  system  Is  type  (111))  and  perform  the  Identical 
analysis,  or  If  the  system  Is  type  (11)  then  there  Is  no  attractive  half  line, 
(see  second  footnote  on  p.  16). 
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